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ABSTRACT 



A method and apparatus are provided that allows initial 
beamforming weights to be determined before timing or 
delay and frequency offset are known. According to one 
aspect of the present invention, the invention includes 
receiving a burst from a terminal at a set of spatial diversity 
antennas, the burst including a known training sequence and 
storing measurements of the received burst from each spatial 
diversity antenna corresponding to receive times within a 
selected analysis window. The invention further includes 
determining weights for each diversity antenna for at least a 
portion of the training sequence of the burst before timing 
for the burst has been determined, applying the determined 
weights to the stored measurements and determining timing 
for the burst based on the stored measurements to which the 
determined weights have been applied. 

31 Claims, 7 Drawing Sheets 
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FIG. 2 
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CALCULATING SPATIAL WEIGHTS IN A aspect of the present invention, the invention includes 

RADIO COMMUNICATIONS SYSTEM receiving a burst from a terminal at a set of spatial diversity 

antennas, the burst including a known training sequence and 
storing measurements of the received burst from each spatial 
BACKGROUND OF THE INVENTION 5 diversity antenna corresponding to receive times within a 

■4 rj« i j c t i j selected analysis window. The invention further includes 

1. meld oi the Invention « . . . . c , . , £ 

determmmg weights for each diversity antenna for at least a 

The present invention applies to the field of spatial portion of ^ training ^encc of the burst before timing 

diversity radio communications systems and, in particular, for me blirst has 5een determined, applying the determined 

to determining spatial weights for a received burst for which 10 weights to thc stored measurements and determining timing 

timing are unknown. for tne burst based on me stored measurements to which the 

2. Description of the Prior Art determined weights have been applied. 
Mobile radio communications systems such as cellular 

voice and data radio systems typically have several base BRIEF DESCRIPTION OF THE DRAWINGS 

stations in different locations available for use by mobile or is ^ present invention is illustrated by way of example, 

fixed user terminals, such as cellular telephones or wireless an d not by way of limitation, in the figures of the accom- 

web devices. Each base station typically is assigned a set of pa nying drawings in which like reference numerals refer to 

frequencies or channels to use for communications with the similar elements and in which* 

user terminate. The channels are different from those of FIG. 1 is a simplified block diagram of a base station on 

neighboring base stations m order to avoid interference 20 which att embodiment of thc invention can be implemented; 

between neighboring base stations. As a result, the user , . , 

terminals can easily distinguish the transmissions received FIG / 2 » a blo f du & am of a remote tcrminal 011 whlch 

from one base station from the signals received from an embodiment of the invention can be implemented; 

another. In addition, each base station can act independently FIG. 3 is a diagram of an exemplary broadcast channel 

in allocating and using the channel resources assigned to it. 25 BCH burst; 

Such radio communications systems typically include a FIG. ^ is a diagram of an exemplary configuration request 

broadcast channel (BCH). The BCH is broadcast to all user CR bu rst; 

terminals whether they are registered on the network or not FIG. 5 is a diagram of an exemplary training sequence in 

and informs the user terminals about the network. In order accordance with the present invention; 

to access the network, a user terminal normally tunes to and 30 FIG. 6 is a flow chart of an implementation of aspects of 

listens to the BCH before accessing the network. It will then the present invention; 

use the information in the BCH to request access to the rin n ■ „ * n _„ . r n , . „ a 

- _ , . „ t M . , - FIG. 7 is a graph oi a search for a least squares error 

network. Such a request typically results m an exchange of re l a tive timing; and 

information about the network using separate control and _„ _ . , _ _ „ 

access channels and ends in the user terminal receiving an 35 ? 15 a S ra P h of a ^Uow-up search for a least squares 

assignment to a particular base station. ^ crror relatlvc Umm %' 

While frequency and timing offset can sometimes be DETAILED DESCRIPTION OF THE 

determined by a user terminal based on the BCH, the initial INVENTION 

request for access is typically transmitted by the user ter- 4Q 

minal without any knowledge of its distance to the base Base Station Structure 

station, its timing delay or the direction to the base station. preserjt i nven tion relates to wireless communication 

In a spatial diversity multiple access system, the base station S y Stems and may be a fixed-access or mobile-access wireless 

can enhance the capacity of the system by determining the network using spatial division multiple access (SDMA) 

position and range to the user terminal as well as any other 45 techn ology in combination with multiple access systems, 

spatial parameters. The timing uncertainty of the arrival time such as thnc division multiple access (TDMA), frequency 

of such request messages are proportional to the round trip division multiple access (EDM A) and code division multiple 

delay encountered by messages traveling between the base access (CDMA)> Multiple access can ^ combined with 

station and the mobile terminal. For systems with a high f reque ncy division duplexing (FDD) or time division 

coverage area per base station, this range and therefore the 5Q duplexing (TDD). FIG. 1 shows an example of a base station 

delay uncertainty may be very large. For example, a range of a wire iess communications system or network suitable for 

of fifteen km results in a roundtrip delay time of around 100 implementing the present invention. The system or network 

microseconds. includes a number of subscriber stations, also referred to as 

In order to accurately resolve the access request and remote terminals or user terminals, such as that shown in 

determine spatial parameters the timing or delay, frequency 55 FIG. 2. The base station may be connected to a wide area 

offset and spatial signature of the user terminal's message network (WAN) through its host DSP 231 for providing any 

are all desired. Typically, it is preferred to accurately deter- required data services and connections external to the imme- 

mine the timing, frequency offset and delay before the diate wireless system. To support spatial diversity, a plurality 

spatial signature is determined. This can create great of antennas 103 is used, for example four antennas, although 

demands on the processing resources of a base station and 60 other numbers of antennas may be selected, 

increase the amount of time required to generate a reply to Aset of spatial multiplexing weights for each subscriber 

the access request. station are applied to the respective modulated signals to 

BRIEF SUMMARY OF THE INVENTION produce spatially multiplexed signals to be transmitted by 

the bank of four antennas. Thc host DSP 231 produces and 

A method and apparatus are provided that allows initial 65 maintains spatial signatures for each subscriber station for 
beamforming weights to be determined before timing or each conventional channel and calculates spatial multiplex- 
delay and frequency offset are known. According to one ing and demultiplexing weights using received signal mea- 



06/23/2004, EAST Version: 1.4.1 



US 6,650,881 Bl 

3 4 

surements. In this manner, the signals from the current active higher processing in the base station, receives downlink 

subscriber stations, some of which may be active on the messages and traffic data from the other parts of the base 

same conventional channel, are separated and interference station, processes the downlink bursts and formats and sends 

and noise suppressed. When communicating from the base the downlink bursts to a transmit controllermodulator, 

station to the subscriber stations, an optimized multi-lobe 5 shown as 237. The host DSP also manages programming of 

antenna radiation pattern tailored to the current active sub- other components of the base station including the transmit 

scnber station connections and interference situation is controllermodulator 237 and the RF timing controller shown 

created. Suitable smart antenna technologies for achieving as 233, 

such a spatially directed beam are described, for example, io ^ nc ( . . . „ . - . At A . n „ 

U.S. Pat Nos 5,828,658, issued Oct. 27, 1998 to Ottersten 10 J* ^J^JTES? ? i T & ™* L f 

et al. and 5,642,353, issued Jun. 24, 1997 to Roy, III et al. 1 T • TT , a t u ^?^tT & ! ^1 

...... . . . , timing signals that are used by both the RF system and the 

The output of the antennas ate connected to a duplexer modem ^ RF controUer 23 3 reads and transmits power 

switch 1«7, which in tins TDD system * a time swnch. Two monitoring and values> conlrols me d , ex / r 107 

possible implementations of switch 107 are as a frequency ^ receives timi parameters and other setlin for each 

duplexer in a frequency division duplex (FDD) system, and 15 burs , from ^ host Dsp 231 

as a time switch in a time division duplex (TDD) system. _ . .... 

When receiving, the antenna outputs are connected via a ^ tnn f l, | : COn i!°Ji er ^ d ^ ator 23 ?> ' eceives . tran ™ lt 

switch 107 to a receiver 205, and are mixed down in analog data from h ° st DSP f ° ur s y mbo3s , at a ^ 

by RF receiver ("RX") modules 205 from the carrier fre- lransmit c ° n [ roller ™* to produce analog IF 

quency to an FM intermediate frequency ("IF"). This signal M ***** ™ h ? h ™ f nt 10 thc ™ transmitter (TX) modules 

then is digitized (sampled) by analog to digital converters 245 ' Specifically, the received data bits are converted into a 

("ADCs") 209. Only the real part of the signal is sampled. ™ m P lcx modulated signal, up-converted to an IF frequency, 

Final down-converting to baseband is carried out digitally. 4 ; times J rT^n^W Y tra f mit ^^S^ 8 

Digital niters can be used to implement the down-converting obtamed from host J^P j?J ? and c ° nverted via digital to 

and the digital filtering, the latter using finite impulse M analog converters ("DACs*) which are part of transmit 

response (FIR) filtering techniques. This is shown as block controllermodulator 237 to analog transmit waveformns. 

213. The invention can be adapted to suit a wide variety of ™ e analo S reforms are sent to the transmit modules 245. 

RF and IF carrier frequencies and bands. The transmit modules 245 up-convert the signals to the 

There are, in the present example, four down-converted transmission frequency and amplify the signals. The ampli- 

outputs from each antenna's digital filter device 213, one per 30 fied transmission signal outputs are sent to antennas 103 via 

receive timcslot. The particular number of timeslots can be ^ duplexertime switch 107. 

varied to suit network needs. While the present example uses j ennma j structure 

four uplink and four downlink timeslots for each TDD * w 

frame, desirable results have also been achieved with three FIG. 2 depicts an example component arrangement in a 

timeslots for the uplink and downlink in each frame. For 35 remote terminal that provides data or voice communication, 

each of the four receive timeslots, the four down-converted The remote terminars antenna 45 is connected to a duplexer 

outputs from the four antennas are fed to a digital signal 46 to permit antenna 45 to be used for both transmission and 

processor (DSP) device 217 (hereinafter "timeslot reception. The antenna can be omni-directional or direc- 

processor") for further processing, including calibration, tional. For optimal performance, the antenna can be made up 

according to one aspect of this invention. Four Motorola 40 °f multiple elements and employ spatial processing as 

DSP56303 DSPs can be used as timeslot processors, one per discussed above for the base station. In an alternate 

receive timeslot. The timeslot processors 217 monitor the embodiment, separate receive and transmit antennas are 

received signal power and estimate the frequency offset and used eliminating the need for the duplexer 46. In another 

time alignment. They also determine smart antenna weights alternate embodiment, where time division diversity is used, 

for each antenna element. These are used in the spatial 45 a transmitreceive (TR) switch can be used instead of a 

diversity multiple access scheme to determine a signal from duplexer as is well-known in the art. The duplexer output 47 

a particular remote user and to demodulate the determined serves as input to a receiver 48. The receiver 48 produces a 

signal. down-converted signal 49 which is the input to a demodu- 

The output of the timeslot processors 217 is demodulated jator 51. A demodulated received sound or voice signal 67 

burst data for each of the four receive timeslots. This data is 50 I s m put to a speaker 66, 

sent to the host DSP processor 231 whose main function is The remote terminal has a corresponding transmit chain in 

to control all elements of the system and interface with the which data or voice to be transmitted is modulated in a 

higher level processing, which is the processing which deals modulator 57. The modulated signal to be transmitted 59, 

with what signals are required for communications in all the output by the modulator 57, is up-converted and amplified 

different control and service communication channels 55 by a transmitter 60, producing a transmitter output signal 61. 

defined in the system's communication protocol. The host The transmitter output 61 is then input to the duplexer 46 for 

DSP 231 can be a Motorola DSP56303. In addition, timeslot transmission by the antenna 45. 

processors send the determined receive weights for each The demodulated received data 52 is supplied to a remote 

user terminal to the host DSP 231. The host DSP 231 terminal central processing unit 68 (CPU) as is received data 

maintains state and timing information, receives uplink burst 60 before demodulation 50. The remote terminal CPU 68 can 

data from the timeslot processors 217, and programs the be implemented with a standard DSP (digital signal 

timeslot processors 217. In addition it decrypts, processor) device such as a Motorola series 5 6300 DSP. This 

descrambles, checks error correcting code, and deconstructs DSP can also perform the functions of the demodulator 51 

bursts of the uplink signals, then formats the uplink signals and the modulator 57. The remote terminal CPU 68 controls 

to be sent for higher level processing in other parts of the 65 the receiver through line 63, the transmitter through line 62, 

base station. With respect to the other parts of the base the demodulator through line 52 and the modulator through 

station it formats service data and traffic data for further line 58. It also communicates with a keyboard 53 through 



06/23/2004, EAST Version: 1.4.1 



US 6,61 

5 

line 54 and a display 56 through line 55. A microphone 64 
and speaker 66 are connected through the modulator 57 and 
the demodulator 51 through lines 65 and 66, respectively for 
a voice communications remote terminal. In another 
embodiment, the microphone and speaker are also in direct 
communication with the CPU to provide voice or data 
communications. 

The remote terminal's voice signal to be transmitted 65 
from the microphone 64 is input to a modulator 57. Traffic 
and control data to be transmitted 58 is supplied by the 
remote terminal's CPU 68. Control data 58 is transmitted to 
base stations during registration, session initiation and ter- 
mination as well as during the session as described in greater 
detail below. 

In an alternate embodiment, the speaker 66, and the 
microphone 64 are replaced or augmented by digital inter- 
faces well-known in the art that allow data to be transmitted 
to and from an external data processing device (for example, 
a computer). In one embodiment, the remote terminal's CPU 
is coupled to a standard digital interface such as a PCMCIA 
interface to an external computer and the display, keyboard, 
microphone and speaker are a part of the external computer. 
The remote terminal's CPU 68 communicates with these 
components through the digital interface and the external 
computer's controller. For data only communications, the 
microphone and speaker can be deleted. For voice only 
communications, the keyboard and display can be deleted. 

Broadcast Channel (BCH) 

In one embodiment, the system of the present invention is 
initiated for each user terminal or remote terminal from the 
broadcast channel BCH which is transmitted as a burst from 
the base station to all potential user terminals. The BCH 
burst, unlike the traffic channel bursts, is transmitted in all 
directions where user terminals may be, typically omnidi- 
rectionally but the specific beam pattern will depend on the 
network. An example of a broadcast burst structure is shown 
in FIG. 3. The BCH communicates enough basic informa- 
tion to enable a subsequent exchange of configuration 
request CR and configuration message CM between the base 
station and the user terminal. The BCH also provides general 
frequency offset and timing update information to all user 
terminals. 

Table 1, below summarizes the content of an example of 
a BCH burst, as shown in FIG. 3. 



TABLE 1 


Duration 


Contents 


10 //sec 


ramp - up 


272 ^scc 


frequency correction training symbols f„ f 2 , . . . , f li6 


256 ftstc 


timing correction training symbols t t> t^ . . . , tj 24 


16 fdscc 


broadcast preamble r lf . . . , r„ 


512 //sec 


information symbols h\, h'a, • . ■ , b'zsc 


10 //sec 


ramp • down 


14 //sec 


inter-burst guard time 



The frequency and timing correction training symbols can 
be set according to any one of many approaches well-known 
in the art. They can also be combined, exchanged with a 
synchronization sequence or eliminated. 

The broadcast information symbols are constructed from 
a 15 -bit broadcast message which is modulated and coded 
into a 256 bit sequence. The number of symbols as well as 
the structure and sequence of transmitted bits can be varied 
to suit a wide variety of applications. The broadcast channel 
information symbols provide the information needed for a 
user terminal to request a configuration message from the 
base station. 
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Each broadcast message is mapped into a broadcast burst 
with the information shown below in Table 2. 



TABLE 2 



5 



10 



Field 


;t Message 

# of Bits 


BStxPwr 


5 


BSCC 


7 


BSload 


3 


Total 


15 



BStxPwr is the effective isotropic radiated power of the 
broadcast message. This number indicates the power trans- 
mitted by the base station taking into account the number of 
amplifiers and diversity antennas available at the base sta- 
tion. 

2Q BSCC is the base station color code, used by the user 
terminal to select training data for uplink bursts and to 
distinguish broadcasts of different base stations. 

BSload is an indication of the amount of unused capacity 
the base station has. 

25 In one embodiment, the network is designed to take 
maximal advantage of spatial division multiple access tech- 
nologies and particularly smart antenna array signal pro- 
cessing. To help maintain reliable spatial channels in an 
extremely dense frequency reuse pattern, the network uses 

30 time division duplex TDMA where uplink and downlink 
transmissions are always on the same frequency. In addition, 
because many user terminals are single antenna and transmit 
and receive omnidirectionally, except for the BCH, an 
uplink burst is always received before a downlink burst 

35 needs to be sent. This allows downlink bursts to be more 
accurately spatially directed. An uplink training sequence is 
embedded in every uplink burst to allow for moderately fast 
frequency hopping despite any decorrelation of the spatial 
channel with frequency. 

40 The frequency hopping sequence may be any one of many 
different sequences well-known in the art. In one 
embodiment, the parameters of the frequency hopping 
scheme are initially unknown to the user terminal. This 
maximizes the flexibility of the network and increases the 

45 flexibility of the user terminal. As explained below, the 
frequency hopping parameters are transmitted to the user in 
the CM burst. 

In one embodiment, the BCH channel is shared by all base 
stations in the wireless communication system. Using the 7 

50 bit BSCC, up to 128 base stations can be accommodated. 
The BCH is part of a time division duplex channel with a 
repeating frame. The channel that includes BCH is a single 
RF carrier frequency used for uplink and downlink. For high 
noise environments or for increased robustness, the BCH 

55 can hop frequencies according to a predetermined scheme or 
be repeated on several different frequencies. The repeating 
frame includes the downlink BCH for each base station, 
labeled BS1 etc. as shown in Table 3 below. The next frame 
includes the uplink Configuration Request CR, labeled CR1 

60 etc. and downlink Configuration Message CM, labeled CM1 
etc. Each frame also includes a number of reserved slots, 
shown as empty boxes below. These slots can be used for 
data traffic, if the broadcast channel is also used for traffic, 
for other control messages or reserved to reduce interference 

65 on other channels in the network. In one embodiment, the 
other traffic channels hop frequencies around and through 
the BCH. The frames are repeated for each respective base 
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station 1 to 128 to build a superframe as discussed in more 
detail below. After the last CM, CM128, the superframe 
repeats and begins again with the next superframe and the 
BCH for base station 1. 

TABLE 3 





Uplink 


Downlink 


Superframe 


Frame 1 








BS1 






1 


Frame 2 


CR1 






CM1 








Frame 3 








BS2 








Frame 4 


CR2 






CM2 
















Frame 255 








BS128 








Frame 256 


CR128 






CM128 






Superframe 
2 


Frame 1 








BS1 






Frame 2 


CR1 






CM1 















10 



15 



20 



Configuration Request Burst Structure 

A configuration request (CR) burst is transmitted by a user 
terminal (user terminal) in order to initiate communications 
or registration with a base station (base station). It is 
transmitted after gathering information about the system by 
listening to (multiple) broadcast channel (BCH) bursts. The 
CR burst, is the first communication from a user terminal to 
the base station and therefore, the user terminal does not 
have any information about its range to its selected base 
station. Accordingly timing, range, and spatial processing 
weights among other things are all unknown to the base 
station. An example of a CR burst is shown in FIG. 4. 

The configuration request burst is composed of several 
fields as shown in FIG. 4 which are listed in Table 4. The 
durations are described in terms of microseconds. In one 
embodiment, a symbol period is 2 microseconds. 

TABLE 4 

Configuration Request (CR) Burst fields 
Duration Content* 



25 



In another embodiment, the BCH is on its own channel 
and CR and CM are on a separate control channel. 
Alternately, one BCH can be provided on a constant fre- 
quency and a secondary BCH can be provided on another 
channel with hopping frequency. The hopping channel is 
described in the CM. 

Registration 

In one embodiment of the present invention, a session 
begins when a user terminal seeks to register with a base 
station. The user terminal does this without knowing its 
relative position to the best station in terms of direction and 
distance (or range) to the base station (BS). Therefore, when 
the user terminal (UT or remote terminal or subscriber 
station) requests registration using a configuration request 
. (CR) burst, the base station receives the transmission from 
the user terminal with a large timing uncertainty and 
unknown spatial parameters or weights when compared to 
traffic bursts that already utilize timing advance directives 
from the base station. 

Before registration, several user terminals may transmit 
CR bursts during the same time slot to register with the same 
base station. CR bursts addressed to different base stations 
may also be received. The base station resolves these 
requests using spatial processing and multiple antennas. By 
combining the multiple antenna measurements 
(beamforming) the base station minimizes the interference 
between signals and maximizes the signal-to-noise ratio 
(SNR) of each received burst. 

Beamforming can be performed in a variety of ways. In 
one embodiment, a training sequence is used with a least- 
squares cost function to determine the weights for a beam- 
former. This allows signals to be classified as desired or 
undesired. The training sequence -based approach uses an 
estimation of the timing and frequency offset of the CR 
burst. This estimation task can be performed by determining 
the least-squares cost for each timing and frequency hypoth- 
esis. This demands great computational resources. Using a 
special sequence design and a beamforming algorithm that 
does not require a search over the entire timing uncertainty 
range alleviates the computational load of the full search. In 
one embodiment, since CR bursts have a high timing 
inaccuracy, a periodic training sequence is dedicated for CR 
bursts alone. 
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The training symbols are allocated 240 microseconds in 
order to allow the signal to be accurately received and 
demodulated before a user terminal has registered and has 
received any knowledge of the system. The training symbols 
are discussed in more detail below. 

The 82 information symbols are constructed from the 
configuration request message using, for example, forward 
error correction coding. In the present embodiment, the CR 
burst is modulated using Jt/2-BPSK modulation in order to 
decrease the peak-to -average ratio of the transmitted wave- 
form. 

The information symbols of the present CR burst are 
mapped out as shown in Table 5, below. Any of the items 
listed below can be deleted and transmitted later during the 
registration cycle or not at all based on the needs of the 
system. CR is scrambled by a function of BSCC ensuring 
that even if there is some interference from CRs sent to 
nearby base stations, the demodulation capture effect of the 
BSCC works out any collisions. In one embodiment, the 
scrambling is performed by taking the encoded bit sequence 
and exclusive OR'ing it with the output of a linear feedback 
shift register. 



TABLE 5 



Configuration Request Message 


Field 


# of Bits 


identity 


8 


utClass 


4 


txPwr 


S 


Tbtal 


17 



identity is a set of unique random bits for each user 
terminal that differentiate simultaneous messages from 
multiple user terminals. Because of the randomness and 
large number of bits, it is unlikely that two user 
terminals will select the same identity code at the same 
time. 
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utClass identifies user terminal capabilities (highest the training sequence. The period of the core sequence is a 

modulation class, frequency hopping capabilities, etc.) compromise between the duration of the search range and 

This sequence identifies the type of user terminal that the possibility of finding sequences of the specified period 

sent the CR. A palmtop digital assistant might have with good autocorrelation and cross correlation properties, 

different capabilities than a desktop computer with a 5 In the present embodiment, constraints on the periodicity, 

fixed dedicated antenna. With utClass, the different mean, autocorrelation and cross correlation are used to 

capabilities can be distinguished. define the core sequences of the training sequence for CR 

txPwr represents the power used by the user terminal to burst acquisition. 

transmit the Configuration Request burst. For example, Periodicity: x(k+P)=x(k), in which x(k) is the periodic 

user terminal power=(2-txPwr-30) dBm. 10 par t of a rc/2-BPSK training sequence and P is the period of 

In the present embodiment, the CR burst includes 106 the sequence. In the present example, P=12 symbols. 

microseconds of extra guard time. As mentioned above, the M ean: the mean of the periodic sequence is constrained to 

CR burst is transmitted by a user terminal without any be zero: 

knowledge of the time required for the burst to travel to the 

base station. In the present embodiment, the desired maxi- 15 p 

mum range from a user terminal to a base station is 15 km. = 0 

The delay of the burst's reception at the base station includes *=J 

the delay in receiving the BCH from the base station. A user 

terminal at range R km receives the BCH burst(s) R/c t , . _ _ , , . . 

seconds after they are transmitted by the base station (c=3x 20 Autocorrelation: The ; out-of phase (symbol spaced) auto- 

10* m/sec). When the user terminal decides to transmit a CR J™ 1 **? * b ° Unded by 1/3 ' WherC T * the phaSe ° feet ° f 
burst, its signal will be delayed by another R/c seconds, so "» out " of -P hase sequence: 
that the delay of the CR burst at the base station relative to 
base station time will be 2R/c seconds. To accommodate 
delays created by a maximum range of 15 km, the guard 25 

period at the end of the CR burst time slot should be at least p " £ 1 / 3, t * 0 

100 microseconds. In addition, in the present embodiment, 
the user terminal randomly delays its CR burst transmission 
from 0 to 9 symbols in order to reduce collisions with other 

user terminals that are sending CR bursts in the same slot 30 Cross correlation: Symbol spaced cross correlation 
and that may be at about the same range from the base between two candidate periodic sequences x(k) and y(k) is 
station. Therefore, a CR burst from a distant user terminal bounded by 1/3: 
may arrive at the base station as much as 118 microseconds 
late as compared to a user terminal that is very close to the 
base station and that does not delay its CR burst transmis- 35 

sion. This accounts for the 106 microseconds of extra guard — 1 <; 15 W 



time that is allotted to the CR burst and the extra margin / £ / £ 

supplied by the 15 microseconds of inter-burst guard time. V *=i \ *«i 
With the inter-burst guard time 3 microseconds of extra 

margin are allotted before the beginning of the next slot. 4 o _ „„_ - 

The present embodiment also includes an additional 15 , At ca f, two ^-BPSK sequences conform to these con- 
microseconds of interburst guard time. This helps to confine stramt f ■ More sequences may be possible if the constraints 
the RF energy from the CR burst to the time-slot allocated arc [ claxed ^ sequences are used as me core 12 
for the CR burst in order not to create adjacent slot inter- fy mbo1 se ^ uenoe 15 re P e ^ d to ^ orm the ^ nodlc 
ference. The interburst guard time provides a total guard 45 Gaining sequence as discussed below. They are: 
time of 121 microseconds to reduce the likelihood that a CR Si»[l J,lj4 -j -1 d ~1 -j -1 -j] 
burst will extend into the next time slot. s 2 =[l,j,l,-j,-l,-j,l,j,-l,-j,-l,j] 

The CR burst is sent on the control carrier, as an example. The bounds on auto and cross correlations help to make 

exactly 2265 /<sec after receipt of a downlink BCH burst. In delayed versions of these sequences to appear partially 

this way, an otherwise uninitialized user terminal can send 50 uncorrelated to an LS beamformer which resolves them. 

CR without any knowledge of the frequency hopping The two periodic sequences can be shared among base 

sequence parameters. As discussed above, the CR burst is stations. In one embodiment, depending on its color code, a 

shorter than a standard uplink time-slot to allow for base station will accept one of the two core training 

unknown tirae-of-flight from the user terminal to the base sequences but not both. The user terminal will be able to 

station and typically arrives late in the uplink time-slot 55 transmit the correct accepted sequence after learning the 

receive window. base station color code from the broadcast (BCH) burst. This 

. . enables a reuse pattern of 2. For reuse patterns greater than 

Training Sequences 2 , some of the constraints on the sequences can be relaxed 

In the present embodiment, known periodic training to yield a larger set of possible periodic sequences to 

sequences are used for the training of base station beam- 60 distribute among the base stations. 

formers in order to extract the desired CR bursts in the In an alternate embodiment, the user terminal selects 

presence of interference. The training sequence includes which of the core sequences it will transmit. This can aid a 

several periodic repetitions of a core sequence followed by base station in resolving two interfering CR bursts. The core 

a marker sequence. Alternatively, the marker sequence may sequence can be selected based on a serial number, product 

precede the repetition of the core sequence. By repeating the 65 number, ID number or other stored number of the user 

core sequence, the search range for a least squares (LS) terminal. For, example if the two sequences Sj and ^ above 

beamformer can be reduced to a single repetition period of are available, the user terminal can read its serial number 
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register. If the serial number is odd it selects sequence si and 
if the serial number is even it selects S 2 . In another 
embodiment, the user terminal can generate a random num- 
ber from 1 to 2. If there are more than two sequences a 
greater range of random numbers can be generated. The 
random number is then mapped into a list of sequences in 
order to select the periodic sequence that is associated with 
the generated random number. 

Transmitting CR Burst and Random Delay 

In one embodiment, the base station transmits broadcast 
channel (BCH) bursts continuously over a pre-established 
known channel. To begin a session with a desired base 
station, an unregistered and unrecognized user terminal in an 
unknown location listens to the BCH bursts. From the BCH, 
it determines the base station with which wants to establish 
communications and selects an appropriate training 
sequence for that base station. It also selects a random 
number D between 0 and 9, and transmits its CR burst after 
delaying it for D symbols. The number for the delay can be 
selected as described above with respect to the periodic 
training sequence. The user terminal can select an internal 
register such as a Serial No. or Product ID No., that contains 
a number greater than 9 and use this as a pointer to select the 
number D. Alternatively, the user terminal can use a random 
number generator to generate a random number between 0 
and 9 and apply this number to select D. 

In one embodiment, the maximum value of random delay 
D is 9 symbols. Therefore, a user terminal can pick one of 
the 10 random delays, {0,1,2, ... ,9} to send its CR burst 
to the base station. In addition, the maximum allowable 
roundtrip delay from the base station to the user terminal is 
assumed to be 50 symbols. Therefore, a CR burst can be 
delayed as much as 59 symbols for a user terminal at the 
maximum possible range. This delay is measured with 
respect to the CR burst of a user terminal that is at zero range 
and that did not delay its transmission. 

The base station receives and resolves the CR burst as 
described below. The roundtrip delay is determined to be 
2(R/c)+D microseconds (since a symbol period is 2 micro- 
seconds where R km is the range of the user terminal. The 
base station, in the Configuration Message (CM) burst, then 
instructs the user terminal to advance its time by x fl =(R/c+ 
D/2) microseconds in its response to the CR burst. The base 
station is unaware of what, if any delay D the user terminal 
has applied to its CR burst. In response to this instruction, 
the user terminal, knowing its random delay D will subtract 
out the random delay from the propagation delay and 
advance its timing by x fl -(D/2) microseconds. 

The random delay of D symbols is particularly valuable in 
a microcellular setting, or when user terminals cluster at 
about the same range from the base station. If the range is 
about the same for two different user terminals, the CR 
bursts from the two user terminals will look identical at the 
base station and may not be resolved. With random delays, 
the CR bursts will be spaced apart making them easier to 
distinguish. 

CR Burst Resolution 

FIG. 5 shows a training sequence of a CR burst for one 
embodiment of the invention. In FIG. 5, x510 is one of the 
two core sequences {s lf s^. x is the vector consisting of 12 
symbols x : to x 12 . In this application, boldface text is used 
to represent vectors and normal text is used to represent 
scalars. The complete CR Burst training sequence is con- 
structed from the core periodic sequence x. Nine repetitions 
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of x; P1,P2, . . . , P9 are followed by a marker sequence M 
505. FIG. 5 also can be viewed from the perspective of base 
station timing. In one embodiment, FIG. 5 shows windows 
that the base station uses to resolve received training 

5 sequences in a CR burst. In this context, the sequences PI to 
M would line up with the windows as shown only if the 
sequences were received from a user terminal that was in the 
same location as the base station, or in other words, if there 
were no propagation delay between the base station and the 

10 user terminal. In this paradigm, a round trip delay margin 
window 520 contains the sequences PI through P5 for a total 
of 60 symbols or 120 microseconds. This accommodates a 
combined delay of 50 symbols due to the round trip propa- 
gation delay between the user terminal and the base station 

15 and 10 symbols of user terminal selected random delay. A 
beamformer analysis window 530 contains core sequences 
P6 through P9 and so is 48 symbols wide. This window will 
capture 4 periods of x even if the CR burst is delayed by up 
to 60 symbols. 

20 Using the observations that fall into the beamformer 
analysis window, a base station can estimate the delay 
encountered by the CR burst modulo -12 symbols since each 
period PI through P9 is identical and 12 symbols long. In 
other words, the measurement in the beamformer analysis 

25 window will be the same for the delays x and t+12, where 
0<t<12. The marker sequence M at the end of the burst can 
be used to resolve this modulo-12 ambiguity. 

In order to demodulate the payload part of the burst, the 
absolute timing information instead of the relative timing or 

30 modulo-12 timing information is desired. As mentioned 
above, the marker sequence (M) at the end of the training 
sequence of FIG. 5 can be used for this task. The marker 
sequence M is set to be M=-x so that one period of the core 
training sequence x augmented with M will be orthogonal to 

35 two concatenated periods of x, i.e., [x, M] [x,x]^=0 in which 
x^Xj, . . . , x 12 ], the core sequence, forms a period of the 
training sequence. 
Since, in the present embodiment, the maximum allowed 

4Q delay is 60 symbols, a roundtrip delay margin window with 
a duration of 5 periods of the core sequence is provided. This 
leaves a duration of 4 periods of the core sequence for the 
beamformer analysis window. Absolute timing is deter- 
mined by determining which of the four possible periods of 

45 the core sequence is the last one received before the marker 
sequence. The beamformer analysis window may contain 
portions of five different core sequences, at least three of 
which are full periods of the core sequence. The first symbol 
of the found sequence is either 12, 24, 36 or 48 symbols from 

5Q the first symbol of the marker sequence. Therefore the 
absolute timing can be determined by testing only five 
hypotheses. In other words, given x, (0^x<12), as the 
modulo-12 delay estimate from the spatial processing 
search, the possible candidates for the absolute delay are: {x, 

S5 x+12, x+24, x+36, x+48}. 

As can be seen from FIG. 5, for any delay less than 60 
symbols, the beamformer analysis window at the base 
station will always contain at least 4 periods of the periodic 
sequence x. As a result, the first symbol, x lf of the first 

60 repetition of the core sequence within the beamforming 
analysis window can be at most 12 symbols away from the 
border of the round trip delay margin window and the 
beamformer analysis window. 
This property of the repeating periodic training sequence 

65 can be taken advantage of when searching for the first 
symbol, x lr of the core sequence x. FIG. 6 shows a flow 
diagram for resolving the CR burst according to one 
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embodiment of the invention. The CR burst is received at the 

antenna array 602 and the measurements from the multiple _ o vJ . ^ , . _ _ _ e 

antennas for the period that covers the beamformer analysis ^baw 
window are stored 604. Then the beamforming weights that 

will yield the minimum least squares (LS) error for each s 

timing hypothesis in the search window of 12 symbols is After computing the covariance matrix R, a matrix L can 

determined. The least squares error results are used to find be found such that 

a coarse relative timing hypothesis 606. Finding this timing Re»L L ff 

includes performing a coarse search for the symbol timing L is termed as the Cholesky factor of R. 

modulo-12 and performing a follow-up search 608. 10 In one embodiment, the location of the first symbol of the 

A weight vector is then determined 610 using the relative ? orc f inin * ™V™* m °/ er f hc firs * \\ s y mbols of the 
... ^f. ... . . j . ,« . A beamformmg search window is searched for using a least- 
timing. This weight vector is applied to the stored measure- b ... , . c ' =* . __ . 

* * L . c u squares processor with a search step of 2/3 symbols. This 

menu to convert the measurements from each antenna ^ ^ lg e of ^ ^ > u Jg 

channel into a suigle channel 612. Then, this single channel is hypotheses and> accordlngly> 18 least-squares (LS) error 
is analyzed to determine the fine timing 614 and the fre- calculations. The vector coarseSearchGrid contains the 
quency offset 616 of the desired signal. The fine timing and delay h yp 0t hesis values for each 2/3 symbol increment (The 
frequency estimates are used to determine a new more ^nit used here is a symbol period): 
accurate weight vector 618 based on the data in the beam- coarseSearchGrid-[0,2/3,4/3,2, . . . , 11 1/3] which con- 
forming analysis window. This improves the quality of the tains me delay values for ±e 18 hypoth eses. 
weight vector and the new weight vector is applied to the For each delay hypothesis k, (1^18), the following 
stored measurements to convert the measurements from computations are performed: 

each antenna channel into a new single channel 620, further Calculate cross-correlation vector p* for each hypothesis: 
improving the quality of the received signal. The processing 

so far cannot perfectly determine the timing information, 25 

because it can only measure time delays modulo-12 />* = 2j r ^ d *^ t+Tk ^ 

symbols, the duration of the repeating core sequence. The ,gT baw 
timing estimate is in the range {0 . . . 12}, i.e., delays of x 

and x+12 symbols will ideally produce the same timing In which r(t) is mc rcC cived signal samples and d(t) is the 

estimate at this point. Using the new weight vector, absolute samp i e s of the desired signal. The desired signal samples are 

timing is determined by finding the marker sequence 622. assumed to be available in a filter bank over-sampled at 24 

This resolves the modulo-12 timing ambiguity. times the symbol rate ^ value Xfr ^ ^ assumed delay for 

When the absolute timing is known, the entire training the k " th hypothesis, i.e., x fr =coarseSearchGrid (k). 

sequence can be identified. Accordingly a new weight vector 35 Back Substitution: The resulting vector p^ applied to 

calculation can be made using the entire training sequence s °l ve for interim vector x* for each hypothesis k, where 

including the marker sequence 624. This redetermined L is tne Cholesky factor as mentioned above: 

weight vector is still more accurate because of the larger Lx^-p* 

number of samples that can be used. It is applied to the Then the least-squares (LS) fit for each hypothesis is 

stored measurements to again convert the measurements 40 computed: 
from each antenna channel into a single channel 626. With 

the new single channel, the CR burst is demodulated and . . ... vl2 vWv 

read 628. = «^«L 

This processing structure yields accurate beamformer 

weights thai cat ,be used to receive the desired signal without 45 Que p6tiodicity of the desiled signal ^ the 

determining perfect Uming for the signal. beamforming analysis window, the cross correlation vector 

computation can be modified to: 

Coarse Search for Relative Timing 

After the CR burst has been received in each of the 50 85 X Mfl + *' + W)M'+24) + r<f+36)1^(r+T l ) 

antennas of the array 602 and the signals have been stored ,6T * 
604. The system can proceed to determine the relative 

timing based on the stored measurements 606. In one m which T ^ a period of the desired signal? { 

embodiment, this process is done using a covariance matrix 55 

computation and a Cholesky decomposition. T^tjcoreSnapPoiitfgtgcoieSnapPoint+ll} 

For measurements that lie in the beamforming analysis Therefore, the above calculation of the cross correlation 

window assume that measurements are collected at 1.5 times vector p k costs one- fourth as many multiplies as it would if 

the symbol rate. In this example, the beamforming analysis the core sequence were not a repeated periodic sequence 

window is of duration 48 symbols or 4 repetitions or periods 60 within the training sequence. Finally, the hypothesis with the 

of the core sequence. Next, assume that the variable core- minimal LS error f^ is picked. An illustration is presented 

Snap Point points to the first snapshot that is in the beam- with respect to FIG. 7. 

former analysis window, and r(t) denotes the snapshot FIG. 7 shows the LS cost function on the vertical axis 702 

(measurement vector) collected at time t seconds after the as a function of the time delay x for the 18 hypotheses on the 

start of the CR burst (relative to base station time). The 65 horizontal axis 704. The horizontal axis is marked in units of 

covariance matrix R can be estimated as follows, where 2/3 of a symbol time in accord with the sampling rate of 1.5 

Tbaw is the time of the beamforming analysis window: times the symbol rate. The selected best coarse search time 
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delay estimate 706 is shown in FIG. 7 at t<=4/3 symbols. 

Note that the delay is measured with respect to the nominal _ v 2 ^ 

time or modulo-12 time at which the CR burst is expected /(To) = ,J£ W M( '* To)l ~*° x ° 

to arrive. This nominal time is named the coreSnapPoint. BAW 
This simphfied computation of the cross correlation vectors 5 

p k based on the periodicity of the training sequence within The beamforming operation results in a single channel 

the beamforming analysis window yields substantially the measurement that shall be named g(t) and defined as fol- 

same accuracy as a regular cross correlation computation. lows: 

Any mismatch of the training sequence in this simplified gCO^w^n^t) 

approach makes only a minimal impact on the performance 10 The actual computations can be simplified because the 

of the beamformers due to the periodicity of the training time range of the beamforming operation does not need to be 

sequence. the whole burst. Note that the frequency estimation (search) 

has not yet been performed. This reduces the search to a 

o ow-up earc single dimension, namely time-delay. The frequency search 

Having made a coarse determination of the relative timing 15 can be avoided because of the duration of the beamforming 

606 a follow-up timing search 608 can be performed. This analysis window, the maximum frequency deviation 

search has one more LS error calculation. The result of this expected in the received signal and the target SNR level at 

calculation can be compared to the best LS error from the the output of the beamformer. When the beamformer analy- 

coarse search. sis window is short enough to make the frequency offset 

The LS error vector f^from the coarse search is treated as 20 thc dcsired si e nal mismatch is not significant 

a circular buffer and concentration is on the midpoint of the enou g h to de gr ade performance. Furthermore, the inability 

range between the best location from the coarse search (that lo ^P 10 ™ thc SNR cvcn b ? ^ an MSE beamformer 

yielded the minimum LS error) and its neighbor that yielded ( that ^ ows the P erf f ct ^variance matrix and steering 

the next best LS error ^ vector) to high levels (above 5 dB) makes the actual beam- 

. , ,', ,„ . , 25 former blind to imperfections in the desired signal (such as 

w fo C °ed rSe dimensional vectors fre<?uency ofifeel) ^ presence Qf frequency offs£t bec£)mes 

were orme . serious at high target SNRs and ignoring large frequency 

updated^ Values-[f 18 ,f j, £^3,^, . . . , f^f^fj] offsets causes the SNR performance to tend to saturate at 

and around 5 dB output SNR. However, this is still not an issue 

updatedCoarseSearchGrid-[-2/3,0,2/3,4/3,3 1/3, 4, . . . , 30 if the CR burst payload can withstand SNRs as low as 0 dB. 

11 1/3, 12] Furthermore, the single channel estimation of frequency and 

Assume that the hypothesis of f 3 (t-4/3) (of the current fine timing that follows will alleviate these problems using 

list) was selected by the coarse search procedure (as shown fewer computations, 

in FIG. 5) because f 3 <f2<f 4 . In the follow-up search . 

P rocedure,theL>SerroriscomputedforT:«l,smcethisforms 35 Hne ^ming bstimator 

the midpoint of the range between the delays corresponding The coarse search and the follow-up search provide a 

to f 2 (t»2/3) and f 3 (x=4/3). This is shown graphically in FIG. timing delay with an effective rate of three times the symbol 

8. rate. The timing resolution can be improved to better 

FIG. 8 shows the LS cost function 802 as a function of the demodulate the payload. Severe intersymbol interference 

18 hypotheses 804 in which the time domain is again 40 can be expected at the worst case timing events (1/6 symbol 

sampled at a rate of 1.5 times the symbol rate. The selected period of timing error). To estimate fine timing 614, the 

best case coarse search time delay estimate 806 is x»4/3 Oerder-Meyr blind timing estimator can be used. This 

symbols. After the additional hypothesis test, the hypothesis estimator is blind in the sense that it does not assume any 

(among the 19 of them) which corresponds to the minimum knowledge of the transmitted symbols. Therefore it is insen- 

least-squares error is again picked. In this example, x 0 =l is 45 sitive to frequency offset errors. The Oerder-Meyr estimator 

the delay value 808 that accomplishes this. can be used as follows: 

Beamforming Weights First > inter P olate g(0 to 3 tirnes the s y mbo1 rate ( from 15 

times the symbol rate) using an interpolation filter, and call 
After the follow-up search, determining the relative tim- 5Q the output g t <t). The duration of the measurements process 
ing delay x to within one third of a symbol period 608, the amounts to 48+12+82=144 symbol periods, this is the guar- 
weight vector w 0 is recomputed to extract the desired signal anteed maximum number of symbols on the right hand side 
as follows: of the coarse timing estimate in the beamforming analysis 

Using the cross-correlation vector p 0 for the best hypoth- window as shown in FIG. 5. Alternatively, the symbols can 
esis: 55 be sampled at a three times rate originally. In such a case, 

only every other sample would be used in the coarse search 

Zand follow-up search, 
leTsAW Next, the interpolated signal is passed through a memo- 

ryless nonlinearity to obtain 

And the corresponding interim vector x 0 : &(0Hfctf)l- 

LXe=p 0 The absolute value nonlinearity generates tones at the 

The weights can be computed by solving the following multiples of the symbol rate. The nonlinearity is approxi- 

back substitution: mated by fitting a polygon. For example, if z is a complex 

L"w 0 »P 0 65 number with representation z-z r +jz„ then |z|-anax(lzj,|z ( j)+ 
With the resulting LS error which is already determined to 0,34xmin (|z r |,|zj). The zero-crossing of the complex sine 

be f(x 0 ): wave at the symbol rate is determined by using a 3 point 
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DFT and this location is termed as a symbol transmission 
instant. Finally, the symbol transmission instant is found that 
is closest to the position identified by the follow-up search 
as the location for the first symbol of the periodic training 
sequence. 

The closed-form estimate obtained using the Oerder-Meyr 
timing estimator is termed as If this estimate is not in the 
range (0,12), it is mapped into this range by adding or 
subtracting 12 to it. 

Coarse Frequency Estimator 

After fine timing is acquired 614, the measurements are 
used in the analysis window together with the desired signal 
to compute the cross correlation function and determine the 
coarse frequency estimate 616 as the peak of the cross 
correlation function. For simplicity, this analysis is done at 
the 1.5 times over sampling rate (instead of the 3 times rate 
used for the fine timing search), however many other sam- 
pling rates can be used. Specifically, for the m-th frequency 
hypothesis, the function freqMatch(m) is computed: 



10 



15 



20 



freqMaichijn) = 



(/■n>,)exp(-)2*/ M /) 



25 



The frequency candidate that maximizes freqMatch(m) is 
selected as the frequency estimate, fcoarsc Frequency can- 
didates (f m 's) are in the set: 

{-2800,-2000,-1200,-400,400, 1200, 2000,2800} 
Alternatively, the frequency uncertainty can be regarded 
as much less and fewer hypothesis tests can be performed. 

More Accurate Beamforming 

Using the frequency and timing estimates l finei f coarse , an 
improved weight vector, w^, can be determined 618 and 
applied to the found core sequences to obtain an estimate of 
the desired burst s(t): 

The more accurate weight vector is obtained by employ- 
ing d^+Tjj^) exp (J2jtf^t) as the desired signal in the LS 
error computations. However, if the corresponding LS error 
is larger than the follow up search based LS error f c , the fine 
timing and frequency estimates are rejected and the esti- 
mates are set to x fine ^z a and f eoorP<r .0, which results in 

Wjw-W,,. 

To further improve the SNR results at high post-copy 
SNR levels, the search can be performed for a fine frequency 
estimate around the coarse frequency estimate at locations 
f COfl r«+[-800,0,800] using an LS fit to multichannel mea- 
surements for t eTg^. Such a search can be included to 
improve the frequency estimation performance for the CR 
burst demodulator. 

Absolute Timing Recovery 

Prior to determining the absolute timing, as opposed to 
relative timing, the more accurate weight vector is applied in 
order to convert the measurements stored from each of the 
antennas into a single channel measurement 620. This 
conversion is still done using the data in the beamforming 
analysis window because only this data is certain to consist 
only of periods of the core sequence. An appropriate pro- 
cedure 622 for determining the absolute delay can be 
implemented as follows. Assume that the modulo- 12 time- 
delay estimate is x, and that g(t) denotes the output of a 
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beamformer (as described above g(t) is a single channel 
measurement) that extracts the desired CR burst. First, the 
candidate starting locations for the concatenated sequence 
[P9,M] (see FIG. 5) for each of the five delay hypotheses is 
determined. For example, if the absolute delay equals x, the 
sequence [P9,M] should start x+36 symbols away from the 
left border of the beamformer analysis window. For the 
hypothesis in which the absolute delay is assumed to be 
equal to x+48, the sequence [P9, M] should start x 84 
symbols away from the left border of the analysis window. 

More specifically, assume that T k denotes the begin time 
for the k-th hypothesis, i.e., T^-x+J6+(12xk), for k-0, 
1,2,3,4. Let v k denote the samples of g(t) that correspond to 
the interval which starts at T k for a duration of 24 symbols 
(two periods). For each hypothesis, compute: 

M*)-hfc-*JT- 

Pick the hypothesis that yields the maximum h(k). 

Assume that hypothesis k 0 is selected. Then, the absolute 
timing estimate is determined as x+(k 0 xl2). 

Due to the orthogonal nature of [x, m] and [x, x], there is 
a five-ary pulse-position modulation (PPM) demodulation 
problem in which the SNR changes depending on the true 
position of the desired vector [x, m]. For example, if the 
actual delay is less than a period, then the payload portion 
of the desired CR burst will appear as noise for the hypoth- 
esis tests 1,2,3,4. On the other hand, if the actual delay is, for 
example 20 symbols, then, the CR burst payload symbols 
will not contribute to the test of a hypothesis in which k=0 
as noise because of the orthogonality of [x, m] and [x, x], but 
they will contribute to the noise component in h(k) for 
k=2,3,4. With this delay condition, the beamforming analy- 
sis window receives the following 48 symbols: 

X5,X 6 ,X7,Xg,Xg,X 1 o 7 xii»X 1 2>Xi,X2^3,X4, . . . X5,Xg,X7,Xg,Xg, 
X 10> X 1 1 >Xl2> X l > X 2» X 3> X 4> 

In this stream, the first and last 12 symbols are written out. 
The spatial processing search will discover the first occur- 
rence of x ( at a delay of 8 symbols. The next step is to 
determine whether the absolute delay is 8,20,32,44,56 sym- 
bols. In order to determine this, the measurements following 
the beamforming analysis window are focused upon. With a 
delay of 20 symbols, the next 40 symbols will be: 

X 5» X 6> X 7> X 8> X 9> X 10> X ll> X 12» X l> X 2> X 37 X 4> X 5> X e> X 7> X 8> X 9> X 10» 
X .11» X 12>~ X 1>~ X 1>"* X 2»"* X 3»~ X 4>~ X 5>~ X 6> — X7,-X8,-Xg,— 
X 10- X lH- X 12> d l> d 2» d 3 ) d 4» d 5.de> d 7» d 8 

in which d, indicates the i-th payload symbol, which is 
unknown. As a result, v 0 =[x, x] and Vj«[x, -x], which 
results in h(0)=0 and h(l)«24 2 . Accordingly, the 
desired signal does not contribute to the decision vari- 
able h(0) when the absolute delay is larger than 12. On 



the other hand, v 2 *=[-x 1 ,-x 2 , . . . ,-x 12 , d A ,d 2 
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is not necessarily equal to zero. The first half of v 2 
equals the marker sequence, which yields the result: 

This differs significantly from: 

h(2) is expected to be higher than h(3) on average because 
of the coherence term (-12) inside the magnitude operator 
and its presence creates a significant competition to the 
correct hypothesis decision value (h(0)) in the presence of 
measurement noise. Therefore, absolute delay errors that are 
equal to +12 should dominate other values of delay estima- 
tion errors. One way to eliminate the self -noise problem is 
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to put the marker sequence in the beginning of the burst and Cholesky factorization performed. The conversion can be 

search for the concatenated sequence [-x, x]. In such a applied to the entire CR burst so that the CR burst can be 

design, self noise will exist only due to the presence of the demodulated as a single channel 628. Alternatively, the 

ramp-up symbols, but this will affect the noise component in absolute timing can be applied to determine the beginning 

only one of the h(k)'s. Other modifications to the ordering s and end of the CR burst and the existing weight vector and 

of the sequences are also possible to suit particular appli- frequency estimates can be used to demodulate the signal. 

cat i ons The choice will depend on the demands of the particular 

As discussed above, the remaining task is to find the svstem and *e condition of the channel. Having determined 

marker sequence. This can be done based on a cross corre- timing and a weight vector, this weight vector can be used 

lation operation with the 1.5 times over-sampled version of 10 10 **** traffic back t0 the remote from which ^ CR burst 

the sequence [x,-x] with v k corresponding to all of the five was received 630. 

hypotheses locations indicating period slips of 0,1,2,3,4. General Matters 

Using the fine timing estimate x^, the time of the final ^ ^ mvcntion ^ eliminate the need to acquire 

beamformer output sample (that arrives at Q which is ^ informatioQ ^ order to perform beam . 

closest to coreSnapPoint+x^ is determined. Let x. denote 15 ^ fc iodkit of t t traini 

the time difference between the chosen measurement at t c ^ ^ tQ ^ tim / delay enc0UD t er ed by 

and the location of actual timing location: ^ desired ^ moduk) _ the period of thc training 

t =r -coreSna Point-x sequence which is less than the total time delay ambiguity 

*"* range. The reception of uplink messages at the base station 

Remember that the beginning of the k-th hypothesis which are not perfectly synchronized is enhanced. 

(k=0,l ,2,3,4) will be at times T^coreSnapPoint+x^. Training based interference cancellation methods can be 

Therefore, the vector v k to be tested for the cross correlation used to enhance the signal quality for such bursts. However, 

for the k-th hypothesis will take the form: the training sequence should be aligned with the incoming 

25 signal, and therefore the resulting search process involves 

v t ={^(' c +i2,x(3+A:)) ) g(f c +i2x(3+*)+2/3)L testing every possible delay with a predetermined resolution. 

*,+12x(3*) + 2x2/3X . . . *« c+1 2x(3^35*2/3) * th ,f dela y f ows ' hypothesis testing pro- 

cess becomes increasingly more expensive with computa- 
which is a vector with 36 entries. tional resources. 

The Vj/s are cross correlated with the desired signal vector 30 Using a periodic training sequence eases the efforts of 
d, which is formed from the knowledge of the training hypothesis testing followed by a timing acquisition 
waveform and the timing difference x,: sequence to determine the perfect timing information. The 

periodic sequence decreases the timing ambiguity to within 
tH^(ooreSiiapPoint+T e +i2>3), j 4i (coreSnapPaiiit+T < +i2x3+ a period of the training sequence, which can be selected 

2/3), . . . , 5XcoreSnapPoint« <!+ l2x3 + 3Sx2y3)] ^ smaU ^ ^ resources are ava ilable to handle the hypothesis 

in which sXt) is the training waveform for the desired signal testing task. For example, if me expected delay is 50 

The magnitude squared cross correlations for each ^)s 9 the f eriod °/ the f* ue » ce 15 10 

hypothesis * computed: ^ mbols j Wlth Resolution of a symbol 10 hypotheses can 

Jr be tested mstead of 50, provided that the penodic training 

h(k)-\\ k <f"| 2 40 sequence can be observed in the search window, regardless 

of the delay. In addition, the search operation, which 

Let k,, denote the index of the largest h(k). Then, the involves cross correlations can be simplified by the number 

absolute timing is determined as: 0 f periods in the search window, if the search window 

contains 5 periods of the training sequence (period with 10), 
^h^ 12 ) 45 the cross correlations can be simplified by five-fold. 
The values are in terms of symbol durations in the above O ncc tne timin S hypothesis is verified, only the number of 
equation with the period of the training sequence being 12 periods of the repeating core sequence that has been skipped 
symbols must Dc determined. To perform this more accurately, the 
Alternatively, symbol spaced sampling can be used in the multiple channel measurements (if available) can be con- 
absolute timing recovery. This can eliminate the mismatches 50 verted to a single channel measurement using the search 
created in the 1.5 times measurements that result in loss of procedure results outlined above. After that a timing acqui- 
orthogonality. With symbol spaced sampling, this operation siti<m sequence that is at the beginning of the burst can be 
will require 5 cross correlations of length 24 and comparing found > and t0 find the number of periods that were missed, 
the squared-magnitudes of the cross correlations. So in the previous example, only periods {0,1,2,3,4} may 
After absolute timing recovery 622, the whole training 55 have been missed in the beamforming analysis window, 
sequence can be used to improve the timing and frequency Accordingly, 5 searches can be performed and these are 
estimates and these estimates can be used to perform one performed on single channel measurements rather than 
more stage of beamforming over the whole training multichannel measurements. 

sequence duration. This can increase the number of symbols The invention eliminates the high cost of search on 
used in spatial processing to 120 symbols from 48 and 60 possibly multichannel data for all possible delays that may 
improve the performance at a cost of computational com- be encountered by the user terminal as a function of range, 
plexity The beamforming would involve first determining a The burst structure and periodic training sequence makes it 
still more accurate weight vector using the full training easier to search for the arrival time of the signal and the 
sequence 626 and then applying this weight vector to periodicity of the training further reduces the cost of corn- 
convert the measurements stored for each antenna channel 65 putations. 

into a single channel measurement 628. In doing so a new In the description above, for the purposes of explanation, 

covariance matrix should be computed and therefore another numerous specific details are set forth in order to provide a 
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thorough understanding of the present invention. It will be determining timing for the burst based on the stored 

apparent, however, to one skilled in the art that the present measurements to which the determined weights have 

invention may be practiced without some of these specific been applied. 

details. In other instances, well-known structures and 2. The method of claim 1 wherein determining weights 

devices are shown in block diagram form. 5 comprises selecting weights that will yield a minimum least 

The present invention includes various steps. The steps of squares error for at least a portion of the training sequence 

the present invention may be performed by hardware of the burst. 

components, such as those shown in FIGS. 1 and 2, or may 3. The method of claim 1 wherein the weights relate to the 

be embodied in machine-executable instructions, which may spatial directionality and interference of the received burst, 

be used to cause a general-purpose or special-purpose pro- 1Q 4. The method of claim 1 further comprising applying the 

cesser or logic circuits programmed with the instructions to determined timing to the burst and redetermining weights 

perform the steps. Alternatively, the steps may be performed that will yield a minimum least squares error for the burst 

by a combination of hardware and software. The steps have using the determined timing. 

been described as being performed by either the base station 5. The method of claim 4 further comprising applying the 

or the user terminal. However, any steps described as being redetermined weights to at least a portion of the burst and 

performed by the base station may be performed by the user 35 redetermining timing using the at least a portion of the burst 

terminal and vice versa. The invention is equally applicable t0 which the redetermined weights have been applied, 

to systems in which terminals communicate with each other 6 - ^ method of claim 5 further comprising determining 

without either one being designated as a base station, a user frequency offset using the at least a portion of the burst to 

terminal, a remote terminal or a subscriber station. wh J c ^. thc rcd L ct «*ined weights have been applied. 

. . , .j , ,20 7. The method of claim 5 further comprising applying the 

The present invention may be provided as a computer redetermined ^ to the burst and ^rtheT d ^ erm f ning 

program product which may include a machine-readable ^ wi]1 ieM a minimum ^ ^ errar for ^ 

medium having stored thereon instructions which may be 5urst ^ lhe redetermined timing, 

used to program a computer (or other electronic devices) to 8 ^ metk)d of daim 5 wherein mmpaiing a C0V ariance 

perform a process according to the present invention. The matrix comprises computing a covariance matrix only for a 

machine-readable medium may include, but is not limited to, portion of the known training sequence, 

floppy diskettes, optical disks, CD-ROMs, and magneto- 9. The method of claim 1 further comprising generating a 

optical disks, ROMs, RAMs, EPROMs, EEPROMs, magnet single representation of the received burst after applying the 

or optical cards, flash memory, or other type of media/ determined weights. 

machine-readable medium suitable for storing electronic 3Q 10. The method of claim 9 wherein determining timing 

instructions. Moreover, the present invention may also be further comprises calculating a cross correlation vector for 

downloaded as a computer program product, wherein the the plurality of different timing hypotheses for application to 

program may be transferred from a remote computer to a the least squares calculation. 

requesting computer by way of data signals embodied in a H. The method of claim 1 wherein determining the 

carrier wave or other propagation medium via a communi- weights comprises computing a covariance matrix against 

cation link (e.g., a modem or network connection). tnc expected signal. 

Importantly, while the present invention has been 12 The method of claim 1 wherein determining t ming 

described in the context of a wireless internet data system for ^^ t t^T^^f fTl^L J fK^f 

4 L1 , , A , , A ., . A c different tuning hypotheses for a portion of the known 

portable handsets, it can be applied to a wide variety of training seauence 

different wireless systems in which data are exchanged 4Q 13 f^^od of claim x whercin dctC rmining timing 

Such systems include voice, video, music, broadcast and comp rises determining timing with respect to a repeated 

other types of data systems without external connections. corc xq}lCDCC 0 f the known training sequence and then 

The present invention can be applied to fixed remote termi- determining timing with respect to the entire received burst, 

nals as well as to low and high mobility terminals. Many of 14. The method of claim 13 further comprising redeter- 

the methods are described in their most basic form but steps 45 mining weights using the entire received burst after deter- 

can be added to or deleted from any of the methods and mining timing with respect to the entire received burst, 

information can be added or subtracted from any of the 15. The method of claim 1, wherein the burst is configured 

described messages without departing from the basic scope in accordance with at least one of a TDMA, a FDMA, a 

of the present invention. It will be apparent to those skilled CDMA and a FDD air interface standard, 

in the art that many further modifications and adaptations 5Q 16. A machine-readable medium having stored thereon 

can be made. The particular embodiments are not provided data representing sequences of instructions which, when 

to limit the invention but to illustrate it. The scope of the executed by a machine, cause the machine to perform 

present invention is not to be determined by the specific operations comprising: 

examples provided above but only by the claims below. receiving a burst from a terminal at a set of spatial 

What is claimed is: 55 diversity antennas, the burst including a known training 

1, A method comprising: sequence; 

receiving a burst from a terminal at a set of spatial storing measurements of the received burst from each 

diversity antennas, the burst including a known training spatial diversity antenna corresponding to receive times 

sequence; within a selected analysis window; 

storing measurements of the received burst from each 60 determining weights for each diversity antenna for at least 

spatial diversity antenna corresponding to receive times a portion of the training sequence of the burst before 

within a selected analysis window; timing for the burst has been determined; 

determining weights for each diversity antenna for at least applying the determined weights to the stored measure- 

a portion of the training sequence of the burst before ments; and 

timing for the burst has been determined; 65 determining timing for the burst based on the stored 

applying the determined weights to the stored measure- measurements to which the determined weights have 

ments; and been applied. 
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17. The medium of claim 16 wherein the instructions the known training sequence and then determining timing 
causing the machine to perform operations comprising with respect to the entire received burst. 

determining weights further comprise instructions for select- 25. The medium of claim 16, wherein the burst is con- 
ing weights that will yield a minimum least squares error for figured in accordance with at least one of a TDMA, a 
at least a portion of the training sequence of the burst. S FDMA, a CDMA, and a FDD air interface standard. 

18. The medium of claim 16 further comprising instruc- 26 - M apparatus comprising: 

tions causing the machine to perform operations comprising a set of spatial diversity antennas at a terminal to receive 

applying the determined timing to the burst and redetermin- a burst, the burst including a known training sequence; 

ing weights that will yield a minimum least squares error for a memory to store measurements of the received burst 

the burst using the determined timing. 10 from each spatial diversity antenna corresponding to 

19. The medium of claim 18 further comprising instruc- receive times within a selected analysis window; 
tions causing the machine to perform operations comprising a processor to determine weights for each diversity 
applying the redetermined weights to at least a portion of the antenna for at least a portion of the training sequence of 
burst and redetermining timing using the at least a portion of the burst before timing for the burst has been deter- 
the burst to which the redetermined weights have been 15 mined to apply the determined weights to the stored 
applied. measurements, and to determine timing for the burst 

20. The medium of claim 18 further comprising instruc- based on the stored measurements to which the deter- 
tions causing the machine to perform operations comprising mined weights have been applied. 

determining frequency offset using the at least a portion of 27. The apparatus of claim 26 wherein the processor 

the burst to which the redetermined weights have been 20 further applies the determined timing to the burst and 

applied. redetermines weights that will yield a minimum least 

21. The medium of claim 18 further comprising instruc- squares error for the burst using the determined timing, 
tions causing the machine to perform operations comprising 28. The apparatus of claim 27 wherein the processor 
applying the redetermined timing to the burst and further further applies the redetermined weights to at least a portion 
determining weights that will yield a minimum least squares 25 of the burst and redetermines timing using the at least a 
error for the burst using the redetermined timing. portion of the burst to whichfihe redetermined weights have 

22. The medium of claim 16 wherein the instructions been applied. 

causing the machine to perform operations comprising 29. The apparatus of claim 27 wherein the processor 

determining the weights further comprise instructions to further determines frequency offset using the at least a 

compute a covariance matrix against the expected signal for 30 portion of the burst to which the redetermined weights have 

a portion of the known training sequence. been applied. 

23. The medium of claim 16 wherein the instructions 30. The apparatus of claim 27 twherein the processor 
causing the machine to perform operations comprising determines timing with respect to a repeated core sequence 
determining timing further comprise instructions for deter- of the known training sequence and then determines timing 
mining a least squares error for a plurality of different timing 35 with respect to the entire received burst. 

hypotheses for a portion of the known training sequence. 31. The apparatus of claim 26, wherein the terminal is 

24. The medium of claim 16 wherein the instructions within at least one of a TDMA, a FDMA, a CDMA and a 
causing the machine to perform operations comprising FDD radio communications system. 

determining timing further comprise instructions for deter- 
mining timing with respect to a repeated core sequence of * * * * * 
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